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Abstract 

Solutions for the fields in a coated cylinder where the core radius is bigger than 
the shell radius are seemingly unphysical, but can be given a physical meaning if one 
transforms to an equivalent problem by unfolding the geometry. In particular the 
unfolded material can act as an impedance matched hyperlens, and as the loss in the 
lens goes to zero finite collections of polarizable line dipoles lying within a critical 
region surrounding the hyperlens are shown to be cloaked having vanishingly small 
dipole moments. This cloaking, which occurs both in the folded geometry and the 
equivalent unfolded one, is due to anomalous resonance, where the collection of 
dipoles generates an anomalously resonant field, which acts back on the dipoles to 
essentially cancel the external fields acting on them. 

Keywords: Folded Geometries, Cloaking, Anomalous Resonance, Superlenses, Hy- 
perlenses 
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1 Introduction 



Analytical solutions have played an important role in understanding the electromagnetic 
response of inclusions to an applied field. In these analytic solutions nothing prevents 
one from substituting seemingly unphysical values of the parameters. For example, for 
a coated spherical inclusion with core radius Vc and shell radius one may substitute 
into the analytic solution for the fields parameter values Tc and Vg with Tc > Ts- Is 
there any physical significance to such solutions? Introducing the novel concept (from 
the viewpoint of classical electromagnetism) of folded geometries and building upon the 



ideas of Leonhardt and Philbin (2006 ) let us first show that "yes there is" . 

Specifically, for simplicity, we analyze in the quasistatic limit the transverse magnetic 
(TM) solution for a coated cylindrical inclusion. In the usual situation, it is filled with an 
isotropic core material having a homogeneous complex dielectric constant Ec and radius 
Tc, embedded in a homogeneous isotropic shell of dielectric constant having radii 
and Ts, with > r^, which itself is embedded in a homogeneous isotropic matrix having 
dielectric constant e^i- The potential V takes values Vc, Vs and V^m in the core, shell, and 
matrix respectively. Each of these are harmonic functions (satisfying AV = 0) within 
their respective domains, except at singularities which we assume are confined to a finite 
set of points in the matrix. At the interfaces they satisfy the boundary conditions 
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(1.1) 

These equations still make mathematical sense if > rg: we look for harmonic poten- 
tials Vc, Vs and V^ defined in the respective regions r < rc, r^ < r < rc and r > r^, 
and satisfying the boundary conditions (11. ip . where now Ec, Eg and garded as 

mathematical parameters entering these boundary conditions. The dielectric tensor field 
£(x) takes values 

£(x) = in the core, 
= —Egl in the shell, 



E^I in the matrix. 



:i.2) 



with the choices of sign here being motivated by the effect of folding of space "back on 
itself, which affects the direction of derivatives. Indeed, flux will be conserved only 
if the radial component of the displacement fleld D(x) = —£{x.)W changes sign, but 
maintains magnitude, at these interfaces: if V ■ D = then one can draw a flow fleld 
for D with arrows and (by conservation of flux) the arrows must reverse direction at the 
interface. The interface conditions (11.11) are compatible with this constraint provided 
£(x) is given by (11.21) . 
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To make physical sense of such a solution we recall the fact that the quasistatic 
equations (and more generally, the equations of electromagnetism) retain their form 
under coordinate transformations. Specifically if V{x.) is a solution to 

V ■ £(x)VV(x) = (1.3) 

and x'(x) is a transformation to a new curvilinear coordinate system, then the potential 
V'{x.') = V^(x(x')), where x(x') is the inverse transformation, satisfies 

V ■ 6'(x') VV'(x') = (1.4) 

where the dielectric tensor, viewed as a contravariant tensor density, has been transformed 
according to the standard formula 

£'(x') = |detA(x)|-iA(x)s(x)A^(x) (1.5) 

in which A = Vx'(x) is the Jacobian, and x = x(x'). The equation (11.41) can be 
reinterpreted as a quasistatic equation in a body with dielectric constant s'(x') in which 
x' = {x[,X2,x'^) are now regarded as Cartesian coordinates. The displacement field and 
the electric field E(x) = — VV^(x) transform to 

D'(x') = |detA(x)riA(x)D(x), E'(x') = [A^(x)]-1E(x). (1.6) 

To turn the unphysical solution in the folded geometry, with rc > Vs, into a physical 
solution we use a coordinate transformation which unfolds the geometry. Consider the 
standard polar coordinates (r, 6) and (r', 6') in the folded and transformed geometries 
respectively. Then the simplest unfolding mapping, as sketched in Fig.l, is given by 
e' = 9 and 

r' = r~^[rs — a^r^ — rs)]r, in the core, 
= Ts — a{r — Ts) in the shell, 
= r in the matrix, 

(1.7) 

where a is a fixed positive constant less than rs/(rc — rs). We emphasize that the pair 
(r, 6) with rc> r > Tg and 27r > 6* > does not suffice to uniquely specify a point in the 
folded geometry: one has to specify whether the point lies in the core, shell, or matrix. 
In a folded geometry it is as if space overlaps itself but without intersection: as one goes 
continuously on a straight line trajectory from the origin, first one moves in the core and 
the radius increases until one encounters the core radius Tc , then one moves into the shell 
and the radius decreases until one reaches the shell radius r<. < rc, where one moves into 
the matrix and the radius increases again. With this definition, the unfolding mapping 
(11.71) is globally a 1 to 1 mapping. 

It is clear from (ll.7p that = rs > = — a{rc — r^). The inverse folding 
transformation x(x') takes the same form as (11.71) with Tc, and a replaced by r^, 
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Figure 1: Sketch of the unfolding transformation ( ]1.7p . where r and r' are the radial 
coordinates in the folded and unfolded geometries. Note that = since the mapping 
is the identity map in the matrix. 



and a ^ respectively, and the roles of r and r' swapped. Using the expression (11.51) and 
the formula for the unfolding map, which in particular implies that in the shell 



x' = -ax + 6x/ Vx • X, A = Vx = (6/r - a)I - 6x (g) x/r , (1.8) 

where 6 = (1 + a)rs, we get expressions for the dielectric tensor in the core, shell and 
matrix in the unfolded geometry 

, _ , _ {b/r - afl + {2ab/r^ - 6Vr^)x O x , _ 

a{b/r — a) 

To be physically reahzable we require that e'^, e'^ and have positive semi-definite imag- 
inary parts, which requires that Ec and have a non-negative imaginary part, while Es 
has a non-positive imaginary part (as can be seen directly from (11.21) and (II. 5p ). In sum- 
mary we see that seemingly paradoxical geometries may be transformed into a physically 
comprehensible form, which may prove an interesting direction for future research. 

When Es = £m the response of the coated cylinder in the folded geometry is equivalent 
to that of a solid cylinder of radius and dielectric constant ^c- The potential in the shell 
in the folded region between r,. and is the same as that in the matrix in this region and 
is the analytic extension of the potential surrounding the solid cylinder provided there 
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are no singularities in this analytic extension- otherwise a solution does not exist. So in 
the unfolded geometry the shell with dielectric tensor s'^lx.') acts to magnify the core by 
a factor of r^/r^ so it responds like a solid cylinder of radius Vc and dielectric constant Ec- 
We call such a shell an impedance-matched hyperlens lens in recognition of the pioneer- 
ing work of Kildishev and Narimanov (2007) who showed that it would magnify fixed 
sources in core, not just in the quasistatic limit, but also for the full Helmholtz equation 
(provided the magnetic permeability was also suitably chosen). Such lenses were first 
considered by Rahm, Schurig, Roberts, Cummer, Smith, and Pendry (2008) as electro- 
magnetic concentrators. Although both groups assumed Ts > Tc, their analysis extends 
directly to the case Vc > Vg. Other hyperlenses with magnifying properties were studied 
by Jacob, Alekseyev, and Narimanov (2006) and Salandrino and Engheta (2006). 

This equivalence is similar to the result of Nicorovici, McPhedran, and Milton (1994) 
who found that a coated dielectric cylinder with radii rg > Vc and moduli Ss = — £m 
would have the same quasistatic response as a solid cylinder of radius r* = r^/vc and 
dielectric constant ^c, i-e. the shell, of dielectric constant = —£m, now known as a 
cylindrical superlens, acts to magnify the core by the factor h = rl/r1. This equivalence 
implied that a line source at radius tq > rg in the matrix would generate a potential 
which appeared like it originated from the line source plus an image line source at the 
radius r^/rg which would he in the matrix when tI/tq > Vg. They found that the actual 
potential in the matrix converged as Eg — to this singular potential at radii greater 
than r^/ro and numerically found that the actual potential developed large oscillations 
at smaller radii. (See, in particular, the sentence beginning with "These fluctuations 
become less pronounced.." above figure 2 in that paper.) To our knowledge this was the 
first discovery of an apparent (ghost) singularity in the field surrounding an inclusion, or 
in effect the first example of perfect imaging (in quasistatics) of a point or line source. 
The regions where the field diverges were later called regions of anomalous resonance 
(Milton, Nicorovici, McPhedran, and Podolskiy 2005). 

In a subsequent development Pendry (2000) made the bold claim that the Veselago 
lens ( Veselago 1967 ) consisting of a slab of thickness d with dielectric constant = — 1 
and magnetic permittivity fig = —1, surrounded by a medium with dielectric constant 

= 1 and magnetic permittivity /im = 1, would behave as a superlens: a line source 
at a distance do in front of the slab, would appear to have an image line source at a 
distance d — do behind the slab. When Eg and fig approached —e^ and — /^m (having a 
very small imaginary part) the actual fields behind the slab converged to these singular 
fields behind the image, but diverged between the image and the slab. There was also a 
seeming paradox (pointed out to GWM by Alexei Efros): if the source was closer than 
a distance d/2 to the lens then the electromagnetic power dissipated in the lens per unit 
time by a constant amplitude source would approach infinity as the loss went to zero. 
This paradox was resolved by Milton, Nicorovici, McPhedran, and Podolskiy (2005) who 
showed that when do < d/2 then the anomalously resonant fields acting on the source 
act as a sort of "optical molasses" against which the source has to do a tremendous 
amount of work to maintain its amplitude. Subsequently it was found that a polarizable 
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dipolar line source or single constant energy line source becomes "cloaked" if it is within 
a distance d/2 oi the slab lens or within a radius \/rfJr\. of a cylindrical superlens (with 
the core having dielectric constant Ec = ^m)- Its dipole moment, and consequently its 
effect on the field outside a certain distance from the lens, becomes vanishingly small. 
The energy generated by a constant energy source, like the energy generated by two op- 



posing sources on opposite sides of a slab lens (Cui, Cheng, Lu, Jiang, and Kong 2005 



IBoardman and Marinov 20061) is effectively trapped within the cloaking region. This 
cloaking was proved (IMilton and Nicorovici 20061) and numerically verified ( Nicorovici, Milton, McPhedra 
to extend to collections of finitely many polarizable dipoles. Also arguments were pre- 
sented (Milton, Nicorovici, and McPhedran 2007) which suggested that a line dipole which 
was "switched on" at time t = in front of a perfect lens with no loss, having Eg = —e^ 
and /is = — /im, would become cloaked in the limit t — oo. On the other hand 



Bruno and Lintner (2007) showed that a dielectric body such as a solid cylinder of fi- 



nite radius in the cloaking region would only be partially and not fully cloaked in the 
limit as the loss goes to zero. One can conclude that a dielectric body is neither perfectly 
cloaked nor perfectly imaged by superlenses (in the limit as the loss goes to zero) if it lies 
within the cloaking region. 

Here we show that anomalous resonance and cloaking extends to folded cylindrical 
geometries, and therefore also to the equivalent unfolded cylindrical geometries. This 



is not too surprising. Leonhardt and Philbin (2006) realized that the solution for the 



electromagnetic fields in the slab superlens can be viewed as the result of an unfold- 
ing of space, and we know that anomalous resonance and cloaking are associated with 
superlenses. 



There are important conceptual differences between the work of Leonhardt and Philbin (2006), 
and our work. In their work the unfolding transformation is applied to empty space, so 
that in the appropriate region one point gets mapped to three points, and a field at 
that point gets mapped to three fields. In this context it is correct, as they do, to take 
transformations of the moduli of the form (11.51) . but without the absolute value around 
the Jacobian of the determinant. In our approach, applied to the idealized "perfect" 
superlens with the full-Maxwell equations, the unfolding transformation is applied to a 
folded geometry, and there is globally a one to one correspondence between points in the 
folded geometry and the unfolded geometry. (The value of x in the folded geometry is 
not necessarily sufficient to specify a point: one also has to specify the manifold on which 
the point lies.) Given empty space one first inserts a fold. In the half of the fold that 
gets mapped to the lens Vx gets replaced by — Vx in the Maxwell equations because 
of the change in handedness of the space and the moduli are negative to ensure that the 
Maxwell equations V ■ D = V ■ B = remain satisfied in any source free region in the 
folded geometry. At a given value of x in the fold the electromagnetic fields take the 
values (E, D, H, B), (E, -D, H, -B), and (E, D, H, B) on the three different manifolds, 
where E, D, H and B are the electromagnetic fields at x in the original empty space. 
Thus the total displacement field density at x is D (and not 3D). When transforming 
the moduli absolute values around the Jacobian of the determinant are needed to ensure 
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that Maxwell's equations remain satisfied in the unfolded "perfect" superlens geometry. 
Our introduction of folded geometries greatly enlarges the class of geometries to which 
one can transform to simplify the analysis of a problem. This simplification is analogous 
to the way one uses conformal transformations to map to a simpler problem. 

For simplicity our analysis [which for the most part only requires minor modifications 
of the analysis of Milton and Nicorovici (2006)] is for two-dimensional quasistatics. Pre- 
sumably analogous results hold for the full (time harmonic) Maxwell equations in three 
dimensional folded spherical geometries, although we have not explored this. Throughout 
the paper we use the symbol = to mean equal by definition, and the symbol ~ to mean 
approximately equal to. 



2 The Green function for a monopole and solutions 
for a dipole in the matrix 

Let us consider the Green function V{x.) for a point source (monopole) located in the 
matrix. Although unphysical (because the net charge associated with the singularity 
oscillates in time) it is mathematically well defined, and useful for deriving the potential 
associated with a dipole. This potential, by definition, takes values Vc, Vs and in the 
core, shell, and matrix which satisfy 

AK = 0, AV; = 0, AKn = -5(x-xo) (2.1) 

in their respective domains, together with the boundary conditions (11. ip . where 5(y) is 
the standard Dirac delta function for a source located at y = 0. The problem of finding 
V{x) can be solved explicitly using power series with respect to the complex coordinate 
z = Xi + ix2, as follows. Note that the Green function for the Laplace equation in is 
given by the formula 



Vo = -^{log{z - zo) + log{z - zo)) = 



p>°^N-E"-'(^)"-E"-'(|^" 

n=l n=l " 



(2.2) 

This is the potential of a point monopole in a homogeneous free space. 

We are looking for a solution K,c,m to the above problem ( (I2.ip and (11. ip ) in the form 
of a power series in each of the three regions: 

oo oo 

n=0 n=0 

oo oo 

n=— oo ri=— oo 

oo oo 

Vr^ = Vo + Y, 4"^)^-" + ^ 5^)^-" (2.3) 



n=l n=l 
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The substitution of these series in the interface conditions (11 .ip yields via the identity 
rd/{dr) = zd/{dz) +zd/{dz) exphcit expressions for the coefficients An'^'^\ Bn'^'^\ 
The formulae for g m can then be found, and are as follows: 



1 
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(2.4) 



where, in accordance with the definitions in Milton, Nicorovici, McPhedran, and Podolskiy (2005) 



we have introduced the real parameters and 6 (not to be confused with the delta func- 
tion) and the complex parameter ?7sc defined via 



6e" 



(£s + £c)(gm + gs) 



£s + £c 



(2.5) 



These expressions for Vc, Vg and Vm are valid both for the cases > and rc> r^. 

In Fig. 2 we show the potential around a monopole when mapped to the unfolded 
geometry. The contrast is evident between the case of a core of dielectric constant 
matching that of the matrix, which is non-resonant in this example, and the case when 

7^ ^m; which exhibits anomalous resonance. Note that in the first case the coated 
inclusion is almost invisible: the equipotentials outside it are nearly circular. 

By letting zq = roe*^° and differentiating (12. 4p with respect to tq and with respect 
to 6o one obtains formula for the potential associated with a dipole at zq oriented in 
the radial direction, and with one oriented in the tangential direction. The potential 
associated with an arbitrarily oriented dipole is of course a linear combination of these 
two potentials and is given by the formulae 



(A;(i) + A;(2))/ro + A;(i)Fe(z, ^o) + k^'^F,{z,z^) 

(A;(i) + fc(2))/ro + k^'^F^^,{z, Zo) + k^'^F^^,{z,z^) + k^'^Fjz, Zo) + k^'^F,^{z,z^) 



ro{l-z/zo) ro{l-z/zo) 



yZ, Zq 



[Z, Zo 



(2.6) 
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1 2 3 4 5 




Figure 2: Numerical computations for the potential associated with a monopole at zo = 6 
in the unfolded geometry (unfolding parameter a = 0.7) with eg = — 1 + 10~^i, {a)6c = 
£in = 1 and (b)£c = 5, £m = 1- In both cases, Tc = 5.4, = 2, rg = fg = 4. 
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where 

= (-F + ik°) /2, = -(F + zfc°)/2 (2.7) 



in which [in accordance with the definition below equation (3.5) in Milton and Nicorovici (2006 )] 

and /c° are the (generally complex) suitably normalized amplitudes of the dipole com- 
ponents which have even and odd symmetry about the line 6 = Oq, and 

Fc{z,Zo) = —, -S{b,hzlzQ) 

ro(es - ec)(em - e^s) 

Zout(^>^0) = -SibMlz^) 

Zm(2:,2;o) = —, — ^S{b,rll{z%)) 

i'm\Z,ZQ} — 1 (^z.sj 



in which 



h 



and the remaining functions are obtained by replacing z and zq with z and zq in (12. Sp . 



These formulae for the potentials agree with the formulae of Milton, Nicorovici, McPhedran, and Podolski 



and (for a dipole not on the xi axis) with the formulae in the supporting online material of 



Nicorovici, Milton, McPhedran, and Botten (2007) (see http://www.physics.usyd.edu.au/cudos/research 



aside from the (irrelevant) additive constant of [k^^^ + /c^^^)/ro. 

It is interesting to see what happens to the potential in the matrix in the limit as 
approaches e^. Specifically, let us suppose that e^, and remain fixed with 

real and positive, and with Ec possibly complex (with non-negative imaginary part) 
but not real and negative, and that £s approaches Sm along a trajectory in the lower half 
of the complex plane in such a way that 5 ^ oo but remains fixed. We set 

V = (2.10) 

When Es is close to £m (12. 5p implies 

es ~ [1 - 2e'"^/{6r])]e^ when ^ e^, 

^ (1 - 2ie-^'^/Vv^)^m when 6^ = 6^, (2.11) 

and so we have 

5 ^ 2£^/{\es- e,n\\r]\) when ^ e^, 

4e^/|es-emr when = ^m- (2.12) 

Thus for large S the trajectory approaches in such a way that the argument of Eg — £m 
is approximately constant. Since the imaginary part of is strictly negative, while the 
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imaginary part of rj is negative or zero, we deduce that is not equal to vr or — vr and 
this ensures that there are no infinite terms in the series fl2.9l) . 

We need an approximation for S{6, w) in the hmit where 6 is very large. From (12. 9p 
we see that when \w\ < h the series expansion for 6S{6, w) converges in the limit 6 ^ oo 
and as a consequence 

^(^,^)^T77 -y (2.13) 

d[n — w) 

When 1 > \w\ > h the terms in the series for S{6,w) first increase exponentially until 
i reaches a transition region where i ^ n in which n is the largest integer such that 
Sh^ > 1 and after this transition region the terms in the series decay exponentially. To 
a good approximation (which becomes better as 5 ^ oo) we have 

£=1 j=—oo 

Since Sh"^ — > 1 as 5 — > cx), upon solving for n in terms of h and S we obtain 

S{S,w) ^ e-'°s-i°g5/i°g^T(«;) (2.15) 

where 

oo j 

T{w) = y -^^7- (2.16) 

j=-oo 

Assuming z is in the matrix, let us first treat the case when Ec 7^ £m. Then as 5 — 00, 
?7sc approaches rj and for h > r'^/\z'zo\, i.e. for \z\ > r^/ro, (I2.13P implies 

2 

lim F,n{z, zo) = Fi„(z, zo) = , (2.17) 

5^oo rQ{zzQ - r^) 



and as a consequence the potential in the matrix, with \z\ > t'^/tq approaches 



1 ^ ^^c 

ro(l-z/zo) ro(zzo-r^) 



1 ryr, 



+ 



ro(l-2;/zo) roizzo-r"^) 



(2.18) 



which, as might be expected, is exactly the same potential which would be associated 
with line dipole outside a solid cylinder of dielectric constant Ec and radius Vc- In the 
unfolded geometry it appears as if shell has the effect of magnifying the core by the 
factor Tc/r'^ = rc/{rs — a{rc — rg)). When the source is located with rg < tq < r'^/vg it 
will look like there is a ghost singularity in the matrix positioned at z = r^/^. When 
r, < \z\ < rl/vQ fl27[5D implies 5S{5,rl/ {zz^)) scales like 5iog(^?/(kl'-o))/(-iogft) ^nd as a 
result this is a region of anomalous resonance with the potential diverging inside it, 
with this same scaling. 
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When = the same argument shows that as, 5 — oo, Fi^i^z, zq) tends to zero 
for \z\ > r^/rg. In fact it converges to zero in a larger region. To see this, note that 
?7sc scales as 1/v^, and as a consequence of which 5T]scS{5,r'^/{z'Z())) scales like 5'^ where 
T = log(rcrs/ro|-2|)/(— log /i), in the region < \z\ < t'^/tq. This converges to zero for 
\z\ > t^/tq, where r# = y^r^, but diverges to infinity (with increasingly rapid spatial 
oscillations) in the region < \z\ < t^/tq. Thus, as 5 — oo, the potential will 
converge for \z\ > r^/ro to the potential associated with a line dipole in free space, while 
diverging to infinity in the anomalously resonant region < \z\ < t'^/tq. 

It is also interesting to consider the limit as Eg approaches —e^n in the folded geometry. 
The results of Nicorovici, McPhedran, and Milton (1994) apply directly to this case, and 
show that the coated cylinder in the folded geometry is equivalent to a solid cylinder of 
dielectric constant Ec of radius r^/rc, which is less than Vg. In particular, in the unfolded 
geometry, the inclusion will be invisible when Ec = e^: presumably such an object acts as 
a lens to shrink the apparent size of any object inside it. One can check that anomalous 
resonance and cloaking do not occur for sources outside the inclusion in this circumstance. 



3 Cloaking of a single polarizable line dipole 

First we present an example which shows that a polarizable line with polarizability a can 
be cloaked when immersed in a TM field surrounding a folded coated cylinder with core 
radius Vc and shell radius < and with cylinder axis xi = X2 = ^- The polarizable 
line is placed along xi = tq and X2 = 0, where tq > rg. Suppose {Ei{xi, X2), E2{xi, X2), 0) 
is the field with the polarizable line absent (but with the coated cylinder present) due to 
fixed sources not varying in the 0:3 direction lying outside the radius when Ec 7^ ^m, 
and the radius r# = -y/r^ when Ec = e^. We assume these sources are not perturbed 
when the polarizable line is introduced. 

Again, let us suppose that Ec and e^ remain fixed and that Eg approaches E^n along 
a trajectory in the lower half of the complex plane in such a way that (5 — >■ 00 but 
remains fixed. Let us drop the field component of the electric field since it is zero for 
TM fields. The field {E^, E2) acting on the polarizable line has two components: 

{E'„E',) = {E, + EIE2 + El), (3.1) 

where 

E, = E,{ro, 0), E2 = E2{ro, 0), El = El{ro, 0), E^ = El{r,, 0), 

{El{x, y), El{x, y)) = (-aMn(a;i, X2)/dx,, -dV,^{xr, X2)/dx2)^ (3.2) 

and Viii{xi,X2) is the (possibly resonant) response potential in the matrix generated by the 
coated cylinder responding to the polarizable line itself (not including the field generated 
by the coated cylinder responding to the other fixed sources). From (12.61) . (12. 7p and (12.81) . 



or alternatively from (2.5), (3.9) and (3.10) of Milton, Nicorovici, McPhedran, and Podolskiy (2005) 

we have 

^n(xi,a;2) = mz) + /r„(^)]/2 + [/°(^) - /°(z)]/(2z), (3.3) 
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where z = xi + 1x2 and for p=e,o 

in which k*^ and k" are the (suitably normahzed) dipole moments of the polarizable hne 
[k*^ gives the amphtude of the dipole component which has even symmetry about the 
xi-axis while k° gives the amplitude of the dipole component which has odd symmetry 
about the Xi-axis ) and in which q = 1 for p=e and q = —1 for p=o. Differentiating fl3.3p 
gives 



El{x,y) = -z[f^^'{z) - f^,:{z)]/2 - [/°'(^) + /°'(^)]/2, (3.5) 



where 



in which 



r^z'^^c r^z^ 



_ dS{5,w) ^ iw^ 
b {d,w) = ; = > ^ . , , „ ■ (3.7) 



These expressions simplify if z is real since then f^'{z) — f^^'{z) = and {El,El) 
{~fin{^)y ~fin'i^))- particular with 2; = ro we obtain 



El \ _ ^^^^ ( k 

'2 



,e 



where 

We will see that \c{6)\ can diverge to infinity as 5 ^ 00, and that when this happens the 
polarizable line becomes cloaked. 

Now if a denotes the polarizability of the line, then we have 

^^^-a(^J). (3.10) 



-k^J --\EI 

This implies 

which when solved for the dipole moment {k'^,—k°) gives 

k" \ f El 
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where 

a, = [a-' ~ ci5)]-\ (3.13) 

is the "effective polarizabihty" . So far no approximation has been made. 

Notice that when \c{6)\ is very large then ~ —l/c{6). So in this hmit the effective 
polarizabihty has a very weak dependence on a. 

To obtain an asymptotic formula for c{S) when 6 is very large we use the asymptotic 
formula (12.131) and (12.151) . Differentiating these gives 

S'iS,w)^— hi\w\<h (3.14) 

d{h — wY 

for \w\ < h while when 1 > \w\> h 

S'{5,w) ^ -[log5/(w^log/i)]e-'°s"'^°s^/^°s''T(^) + e-^°s'"^°s<5/iog/^2-/^y^) 

^ -[Iog5/(w^log/i)]e-^°s«'iog5/iog?»2-^^)^ (315) 

where T'{w) = dT{w)/dw and in making the last approximation in (13.151) we have as- 
sumed that I log5| is very large. Let us first treat the case where is fixed and not equal 
to Era and To < Tc- Then we have r/sc ~ V and substituting these approximations in (13. 4p 
and (13. 9p and keeping only the terms which are dominant because S is very large gives, 

for r^/ro > \z\ > Tg, 

which implies 

, ^ -gfcPr7e-'^log^ ^_i„g(,g/(,,^)i,g^/,,g,^ , 
ztq log h 

and 

^ !^^ll^g-21og{../.o)logVlog/.2^( 2/ 2)_ (3_ig) 
Tq log Al 

We see that \c{S)\ — > cxd as 5 — cx) when < Vc- Thus for a polarizable line dipole inside 
the radius the "effective polarizabihty" approaches zero in the limit S ^ oo. When 6 
is very large from (I3.12p and (I3.13P we have 

F ^ -Ei/c{6), k° ^ E2/c{6). (3.19) 

Thus for z in the annulus r^/rg > \z\ > the potential associated with the polarizable 
line has, from (13.161) . 

/f„ ^ Ei5'°sW^°)/'°sVoT(r,V(ro^)) \ogh/{T{rl/rl) log 5). (3.20) 

Similarly in this annulus we have 

/° ^ E25'°*^(^/^'°)/'°sVoT(r2/(ro^)) \ogh/{T{rl/Tl) log 5). (3.21) 
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For z outside the radius r^/ro the potential due to the polarizable hue dipole is approx- 
imately given by (12. 181) and converges to zero because h!^^^ and /c*^^-' vanish as 5 — *• cxo. 
We avoid the technical question of what happens when \z\ = r^/rg but presumably the 
potential also converges to zero there. 

Thus as 5 — * oo the potential in the matrix due to the polarizable line dipole converges 
to zero in the region r > tq but diverges to infinity with increasingly rapid angular 
oscillations for Vg < r < tq. (This is to be contrasted with the potential in the matrix 
associated with a line dipole having fixed k'^ and k°, which as can be seen from fl3.16p 
diverges to infinity in the much larger region < r < r'^/rQ.) A simple calculation shows 
that in the shell the potential associated with the polarizable line similarly converges to 
zero for r > tq but diverges to infinity for < r < ro, while in the core the potential 
associated with the polarizable line converges to zero everywhere. 

It is instructive to see what happens to the local field E2) acting on the polarizable 
line as 5 ^ 00. From ([31]), (l332|) and (IXTSD we see that 

E'^.E.,ciS)k'.E.^-Sm_.^^ (3.22) 

goes to zero as 5 — 00, and similarly so too does This explains why the "effective 
polarizability" vanishes as 5 — > 00: the effect of the resonant field is to cancel the field 
{E^, E2) acting on the polarizable line. 

Suppose the source outside is a line dipole with a fixed source term {kf, k°) = (kf, 0) 
located at the point (ri, 0), where ri > > > r^. When ri is chosen with r^/rg > ri > 
Tc the polarizable line will be located within the resonant region generated by the line 
source outside. One might at first think that a polarizable line placed within the resonant 
region would have a huge response because of the enormous fields there. However, we 
will see that the opposite is true: the dipole moment of the polarizable line still goes to 
zero as 5 00. From (13. 8p . (13. 5p and (13.170 . with tq replaced by ri, the field at the 
point (ro, 0) when the polarizable line is absent will be 



El = ci{S)kl, E2 = 0, (3.23) 

(3.24) 



where 

C,(5) ^ ^g"^ ^ g- log(rg/(ro n)) log S/ log / ^^^^^^y 



roTi log h 

This and (I3.19P implies the polarizable line has a dipole moment 

k' ^ -EJc(5) ^ -cA5)kUc(5) ^ _ ^0^(^s/(^O^l)) ^log(n/ro)/logfe^e .3 25) 

riT{r^/r^) 

So k^ scales as S^"s(^"^/'^o)/iogh ^j-^jcfi goes to zero (since h < 1) as 6 ^ oo but fairly slowly 
when ri and ro are almost equal, i.e. both close to Tc. 

If the source is outside the critical radius rcrit = t^/ts then there are no resonant 
regions associated with it and k^ will scale like l/c(5), i.e. as (52iog{r-c/r-o)/iogh^ j^g^ which 
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goes to zero at a faster rate as 5 — oo, but still slowly when ro is close to Tc. On the other 
hand when tq is close to we have r^/ro ~ 1/Vh and this latter scaling is approximately 
6'^/log6 ~ —e'^/loge", where e'^ is the imaginary part of ^s, which is quite fast. 

The asymptotic analysis is basically similar when Ec = and tq < = y^rvr^. Then 
T]sc ~ -ie-'^/^/V6 and from ([32]), ([33]), (I2A5]) . and (IXTSD we have for r^/ro > |^| > Ts 
that 

^ igFe^'^/2e[^°§^-'°s('^^'^=/"o)]^°*5^/'°*5V-iT(r3V(ro2)), (3.26) 

and 

Tq log 

When all the sources lie outside the critical radius Tc so they do not generate any resonant 
regions in the absence of the polarizable line, both and k° will scale as l/c{6), i.e. as 
^iog(rcrs/rg)/iog/i^ jQg^^ as 5 CX3. When ro is close to Vs we have rj-s/rQ ^ l/\/h and 
this latter scaling is approximately 1/ {\f5\og5) ~ s'l/ loge" which is the same as when 
7^ ^m- By substituting 03.191) in (13.261) we obtain 

/fjz) ^ -Ei?e''^/2g[log.-log(rcr./ro)]log5/logV-lT(r2/(roZ))/c(5) 

^ ^^5i°g(-/n))/i°gVoT(r2/(ro^)) log /i/(T(r3Vr2) log 5), 

(3.28) 

which coincides with (I3.20p . Likewise (I3.21|) still holds. By similar arguments applied 
to Vc and it follows that as 5 ^ oo the potential V diverges with increasingly rapid 
oscillations in the core in the region Tc > r > rcrg/ro, in the shell in the two regions 
Tg < r < ro and r,, > r > r^rs/ro, and in the matrix in the region Ts < r < ro. Outside 
these regions it converges to the potential generated by the fixed sources. 

It is possible to get any cloaking radius between and Tc if we let Sc depend on 5, 
so that — £c scales as and Sm — Ss scales as 5"^^^ , where /5 is a fixed constant 
between and 1. Then ?7sc will scale as and c{5) will scale as 5"^ log 5 with r = 
log(r^~^'^rg^/rQ)/(— log/z) and so the cloaking radius will be rl~^r^. Since [based on the 



results of Milton, Nicorovici, McPhedran, and Podolskiy (2005) and Bruno and Lintner (2007)] 



dielectric bodies located in the cloaking region are not perfectly imaged, it is not suffi- 
cient that Ec, Ss, and £m be arbitrarily close to each other to ensure perfect imaging of 
a dielectric body which lies inside the radius Vc- Similarly, for the standard cylindrical 
quasistatic superlens, it is not sufficient that —Eg, and £ni be arbitrarily close to each 
other to ensure perfect quasistatic imaging of a dielectric body which lies inside the ra- 
dius r* = rl/Tf.. Also a slab lens of thickness d and permittivity Eg separating two media 
with permittivities Sm and Sc will not necessarily provide a good quasistatic image of a 
dielectric body which lies within a distance d of the slab, even when 8^ —Eg, and em are 
arbitrarily close to each other 
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4 A proof of cloaking for an arbitrary number of 
polarizable line dipoles 

The concept of "effective polarizability" does not have much use when two or more 
polarizable lines are positioned in the cloaking region since each polarizable line will also 
interact with the resonant regions generated by the other polarizable lines and if the 
polarizable lines are not all on a plane containing the coated cylinder axis then these 
interactions will oscillate as 5 —>■ oo. However we will see here that nevertheless the 
dipole moment of each polarizable line in the cloaking region must go to zero as 5 ^ oo 
and in such a way that no resonant field extends outside the cloaking region. This is not 
too surprising. Based on the results for a single dipole line we expect that a resonant 
field extending outside the cloaking region would cost infinite energy, and the only way 
to avoid this is for the dipole moment of each polarizable line in the cloaking region to 
go to zero as 5 ^ oo. 

Here we limit our attention to the cylindrical lens with the core having approximately 
the same permittivity as the matrix. Also to simplify the analysis we assume the core 
(but not the matrix) has some small loss. Specifically we assume 

£m = 1, es = l-in, ec = l+i'jK, (4.1) 

with K and jk having positive real parts and approaching zero in such a way that 7, 
which could be complex, remains fixed and (j) given by (12.51) also remains fixed. In this 
limit (12. 5p implies [k + jk)!^ ^ 4/(5e*'^) and since k and 'jk have positive real parts we 
deduce that is not equal to tt or — tt. Solving for n we see that 



K ^ 2e-''^/VV5(l + 7) (4.2) 
The potential in the core due to a single dipole in the matrix at zq is given by ( 12. 6p and 

(ED. 

If there are m dipoles at zi,Z2,---,Zm (where Zi 7^ zj for all i 7^ j) all in the matrix 
then, by the superposition principle, the potential in the core is 



where for £ 7^ 



K = '^(^Pz' + Bf'z'), (4.3) 



= TTE% tc^'/'-M- = EC^f (4.4) 



in which rj = \zj\ and 



d[es - ecj{l - Es) 
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depends on 6 through the dependence of Ss and Ec on 5 but tends to e**^ as 5 — oo. 

Let us suppose the dipoles positioned in the matrix at zi,Z2,---,Zg with I < g <m are 
in the cloaking region, while the remainder of the dipoles are outside the cloaking region, 
i.e. 

< r# for all j < g, > ^# for all j > g, (4.6) 

where we allow for the special case where some of the dipoles have \zj\ — r#: as we will sec, 
these are also cloaked. We do not specify how the set of dipole moments {ki, k2, km} 
depends on S except that: 

• We assume that each dipole outside the cloaking region has moments which con- 
verge to fixed limits as 5 ^ oo 

lim {kf\S),kf\5)) = ikf^,kf^) for all j > g. (4.7) 

The dipole moments kf\5) and kf'\5) inside or outside the cloaking region are 
assumed to depend linearly on the field acting upon them, since non-linearities 
would generate higher order frequency harmonics. Some of them could be energy 
sinks, although at least one of them should be an energy source. 

• We assume that in the unfolded geometry the energy absorbed per unit time per 
unit length of the coated cylinder remains bounded as 5 ^ oo, as, for example, 
must be the case if the line sources only supply a finite amount of energy per unit 
time per unit length. We let VFmax be the maximum amount of energy available per 
unit time per unit length. It is supposed that the quasistatic limit is being taken 
not by letting the frequency uj tend to zero, but instead by fixing the frequency uj 
and reducing the spatial size of the system and using a coordinate system which is 
appropriately rescaled. 

We need to show that, because the energy absorption in the core remains bounded, 
the dipole moments in the cloaking region go to zero as 5 ^ oo and the resonant field 
does not extend outside the cloaking region, r < This is certainly true when only one 
polarizable line is present but as cancellation effects can occur (the energy absorption 
associated with two line dipoles can be less than the absorption associated with either 
line dipole acting separately) a proof is needed. 

To do this we bound k^^^ and k^'^^ for any given i < g using the fact that the energy 
loss within the lens is bounded by W^max- If M^c = Wc{S) represents the energy dissipated 
in the core in the unfolded geometry, then we have the inequality 

Wc = {uj/2) / r'dr' / d9'E'{x') ■ Imag(£')E'(x') 
Jo Jo 

= {u/2)el I rdr I d^E(x) •E(x) > {u/2)el I rdr / dO Ei{z)Ei{z)i4.8) 

Jo Jo Jo Jo 
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in which Imag denotes the imaginary part, e'^ = Imag(£:c) and Ei{z) is the xi component 
of the electric field in the core in the folded geometry given by 



(4.9) 



where the derivative dVc/dxi is calculated by substituting z = xi + 1x2 in (14. 3p . 

Substituting this expression for the electric field back in (14 .Sp and using the orthog- 
onality properties of Fourier modes we then have 



m—l 



e=n-m+l k=0 

where the last identity is obtained using fl4.4p with the definitions 



'j) > 



(4.11) 



in which k = n — I, n>m + l remains to be chosen, and i < g. From (14. lip it follows 
that U = Mu and V = Mv, where M is the Vandermonde matrix 



M 



/I 1 1 

zi/ri Z2/ri zs/ri 

{zi/vif {z2/rif {zs/vif 



V(^i/r,)— 1 {z2/r,r-' (zs/r,) 



1 \ 

{Zm/Tif 

[zJr,T-^j 



(4.12) 



From the well known formula for the determinant of a Vandermonde matrix it follows 
that M is non-singular. Therefore there exists a constant Q > (which is the reciprocal 
of the norm of and which only depends on i, m and the Zj) such that |U| > Cj|u| 
and |V| > Cj|v|, implying 



|Up + |Vp > aXH' + |vp) > c\{\u.\' + = c^(|A:«P + \krY)/rl (4.13) 
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Next we need to select n and find a lower bound on bk which is independent of k. Let 
s = — log 5/ logh (so Sh^ = 1) and take n as the largest integer smaller than or equal to 
sson + l>s>n. Then since r^, < rg < and = r^/i^/^ we have 

(r,/ri)2"52 > (^^/^.)2.^2 _ ^2^-2 log(r,/r,)/ logh, _ ^1/2^-2 log(r#/n)/ log h_ (^^^^^^ 

Also the following inequalities hold for m — 1 > k > 

I = 5h' < 5K^ < dh''-^ and dh''-^ < 5h'-^-^ < dh'-"^ = h-"^ . (4.15) 
So it follows that 

\l + 5h''-^e'^\<a= max ll + te^'^l, (4.16) 
and a is independent of 5. From the bounds 04.141) and 04.161) we deduce that 

hk > (s-m)(r,/r,)2^5i/2^-2i°s(^'#A'')/i°s'^|V'(5)|V«' 

> -[(log5/log/i)+m](5i/2^-2i°s(^'#/'-')/i°s'^|V;((5)|Val (4.17) 

Combining inequalities gives 

2Wjuj > <1^('^)l'^/^ (iogj^^iog;,)^-2iog(r.^/n)/i°gfe(|U|2 + |Vp) 
a''(— log h) 

> ""'l^t^^^^'^'f {logs + m\og (4.18) 
a''rf[ — log n) 

in which the real positive prefactor has the property that 



p, ^ hm " = ^ ; Real(e-^/^7/Vl + 7) (4.19) 

<5^oo a^r|'(— log n.) a'^rf{—logh) 

is strictly positive, where Iieal{w) denotes the real part of w. So there exists a 6q such 
that, for all S > 6q and all z < (?, 

4!|fpM>p/2. where p. mm p. >0, (4.20) 

a"^rj^(— log Ai) «<9 

and such that ^ 

log 5 + m log /i > - log 5 > -2 log /i, (4.21) 

which, in particular, ensures that n > m + 1. So we conclude that 

+ < 2(5'°s(''#/"')/^°s^V2Wc/(^plog5). (4.22) 

which, since log (r^/r^)/ log /i is negative, forces the dipole amplitudes k^^^ and k^^^ to go 
to zero as 5 ^ oo (even when = r^) because Wc = Wc{S) < PVmax. 
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Now the superposition principle implies that the potential at any point z in the matrix 

is 

m 

V(^z) = J2 kf^vl'\z) + kfv}'\z), (4.23) 

where vj^\z) (or vj'^\z)) is the potential in the matrix due to an isolated line dipole in 

the matrix at the point Zj with k^^^ = 1, /cj^"* = (respectively with k^^^ = 0, fcj^^ = 1). 
Now according to the analysis at the end of section 2 (which is easily extended to the 
case treated here where Ec depends on 6 as implied by (14.11) and (14. 2p ) it follows that for 
z in the matrix with \z\ > maxjrg, r^/r^}, 

Inn V}'\z) = Vj'\z) ^ n \ ^ 

<5^oo rj{l~z/Zj) 

\hn V}'\z) = Vj'\z) = . . (4.24) 

Also, as shown in the analysis at the end of section 2, if r'^/rj > \z\ > rg, then vj^\z) 

and vj^\z) diverge as where r = \og{rcrs/rj\z\) / {— log h) . If Zj is outside the cloaking 
region (i.e. j > g) then r'^/rj will be less than r#. So using the well known fact that 

lim e{5)f{5) = eo/o, where Cq = lim e{5), fo = lim /(5), (4.25) 

O^OO 0— >0O 0— >oo 

it follows that for all \z\ > r# and all j > g 

^V^{z) = k^V^[z) 

J1)t/(2)/^A _ J1)t7(2) 
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lim k'I^Vl'>{z) = k'I^Vl 

lim k';^V^^\z) = kf^V^'\z) (4.26) 



+ 00 



If Zi is inside the cloaking region (i.e. i < g) and \z\ > r'^/vi then (I4.24p . (14.251) and 
the fact that and tend to zero implies that kl^^v}^\z) and k^'^^V^'^\z) will 

tend to zero. For r'^/vi > \z\ > we have that v}'^\z) and v}'^\z) scale as 6'^ with 
T = \og{rcrs/ {ri\z\)) /{— log h) while from (I4.22p k^^^ and k^'^^ scale at worst as 6~^/{\og6) 
with t = log(r#/rj)/(— log/i). So their product kl^'^v}^\z) or A;j-^^V^*'^''(2)will scale at 
worst as 6'^~'^/log6 where t — t = log(r#/|2;|)/(— log/i). This goes to zero as 5 — oo 
when \z\ > By taking the limit 5 ^ oo of both sides of (14.230 we conclude that 

m 

}ij^V{z) = J2 [kfjv}^\z) + kf(^y}^\z)] for all \z\ > r#, (4.27) 
i=9+i 

which proves that the coated cylinder and all the line dipoles inside the cloaking region 
are invisible outside the cloaking region in this limit. 
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In this proof we have assumed that the dipole positions Zj are independent of 6. If 
they depend on 5 and \zi{5) — Zj{5)\ is not bounded below by a positive constant for 
all i ^ j then it is an open question as to whether cloaking persists. At least in some 



cases it may persist since Nicorovici, Milton, McPhedran, and Botten (2007) show that 



"polarizable" quadrupoles are cloaked. 

5 Numerical examples of cloaking of collections of 
polarizable line dipoles 

Due to the mathematical equivalence between the analysis for the coated cylinder in 
the cases > Ts and rc < Vg, we can use the same numerical tools here as were em- 



ployed in the paper (Nicorovici, Milton, McPhedran, and Botten 2007) to solve for the 



fields in the folded geometry. Then we use the unfolding transformation (11.71) to obtain 
results for the potential in the unfolded geometry, where the permittivity in the shell 
is anisotropic (with a positive definite imaginary part) and given by (II. 9p . We have 
prepared three animations illustrating the cloaking action, one for a pair of polarizable 
dipoles in a uniform external field, and two others for a set of six polarizable dipoles 
arranged on the vertices of a hexagon. (These animations show snapshots of the po- 
tential distribution in space, for a sequence of equilibrium solutions, as discussed by 



Nicorovici, Milton, McPhedran, and Botten (2007)). We present here in Figs. 3 and 4 



images from each animation. 

Fig. 3 shows the potential associated with two polarizable dipoles, of which one is 
inside the cloaking radius and the other outside it. The resonant region touches the 
cloaked line dipole, and quenches the field acting on it. As in the previous study 



(Nicorovici, Milton, McPhedran, and Botten 2007), the resonance develops first on the 



shell-core boundary, before developing on the shell- matrix boundary (see movie 1). 

Fig. 4 shows two frames from movies 2 and 3, and compares the cloaking of a set 
of six polarizable dipoles for two values of the imaginary part of ^s- As can be seen 
from the first figure, an imaginary part of 10~^ is not sufficient to ensure cloaking of 
the two dipoles closest to However, good cloaking of all six dipoles is achieved 



for an imaginary part of 10 As in the papers of Bruno and Lintner (2007) and 



Nicorovici, Milton, McPhedran, and Botten (2007), it appears that cloaking becomes more 



difficult as the number of polarizable particles in the collection increases, and becomes 
more effective as the particles move more deeply into the cloaking region. 
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Figure 3: Numerical computations for the potential associated with a pair of polarizable 
dipoles (polarizability a = 2 located at points zi = 2.816 — 4.3282 and Z2 = 3.755 + 4. 828z 
in the unfolded geometry (unfolding parameter a = 2) with Es = + Ec = e^a = 1, 

Tc = 8, = 2, Ts = Tg = 4. The dashed hue denotes the cloaking radius, at r# = y/VcTs — 
5.657. Note that one dipole is outside the cloaking region, while the other is inside. 
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